Fermionic superfluid properties in a one-dimensional optical lattice 
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We discuss various superfluid properties of a two-component Fermi system in the presence of a 
tight one-dimensional periodic potential in a three-dimensional system. We use a zero temperature 
mean field theory and derive analytical expressions for the Josephson current, the sound velocity 
and the center of mass oscillations in the BCS-Bose Einstein condensation crossover region. 



I. INTRODUCTION 

Recent experiments with trapped alkaline Fermi gases 
have shown promising directions toward understanding 
the many body physics in various disciplines. These 
systems have became a focuss of attention due to their 
easy controllability. For example, the controllability of 
interaction strength using a magnetically tuned Fesh- 
bach resonance [l| and of the spatial structure using 
counter propagating laser beams allow one to study a 
wide range of quantum mechanical phenomena. One 
such phenomena is a dissipationless flow at low temper- 
atures. The focus of this paper is the analytical study of 
dissipationless transport properties of a stack of weakly- 
coupled two-dimensional (2D) Fermi systems. Such a 
system can be created by applying a relatively strong 
one-dimensional (ID) optical lattice potential to an ordi- 
nary three-dimensional system. One natural analogy of 
this system is the weakly coupled copper-oxide planes of 
high temperature superconductors where Bose condensed 
paired electrons flow dissipationlessly in the crystal lat- 
tice. 

In this paper, we derive analytical expressions for the 
energy and the density in the limit of weakly coupled 
layers using a zero temperature mean field theory. Us- 
ing the energy and the density expressions together with 
gap equation, we then derive analytical expressions for 
various dynamical quantities in the Bardeen, Cooper 
and Schrieffer-Bose Einstein condensation (BCS-BEC) 
crossover region. A similar discussion of tunneling dy- 
namics using a path integral approach can be found in 
ref. [2|- A hydrodynamic approach in the weakly inter- 
acting BCS limit and a Bogoliubov-de Gennes approach 
in the unitarity limit can be found in ref. i3j and ref. [J] 
respectively. First, we calculate the superfluid oscillation 
(Bloch oscillation) amplitudes due to the linear force on 
atoms in the presence of Bloch bands and find that the 
amplitudes are larger for weak periodic potentials and for 
weak attractive interactions. Second, we calculate the 
sound velocity along the direction of ID lattice potential 
and find that the sound velocity strongly depends on the 
lattice height but weakly depends on the interaction for 
tight periodic potentials. Finally, we calculate the dipole 
oscillation frequencies in the BCS-BEC crossover region 
in the presence of a harmonic external trapping potential. 



In the absence of a periodic potential, dipole oscillations 
are undamped and the oscillation frequencies are inde- 
pendent of the two-body interaction (generalized Kohn's 
theorem |5|). However, in the presence of a ID periodic 
potential, we find that the dipole oscillations are damped 
and strongly depend on the two-body interactions. Such 
damping for a ID Bose gas and a non-interacting Fermi 
gas have already been observed loj. 

The paper is organized as follows. In section II, we 
present the theoretical formulation. In section III, we 
present the results and discussion followed by the con- 
clusion in section IV. 



II. FORMALISM 

We consider an interacting two-component Fermi 
atomic gas trapped in one dimensional optical lat- 
tice created by standing laser waves. The ID opti- 
cal potential modulated along z-axis has the form V = 
sEjisin^ {2Trz / X) . Here A is the wavelength of the of the 
laser beam, Er = Ti? (2it / \Y / {2M) is the recoil energy 
and s is a dimensionless parameter which accounts the 
intensity of the laser beam. The periodicity of the ID 
optical lattice along z-direction is d = A/2. When the 
laser power is large, the atomic system forms a stack of 
weakly coupled 2D planes. We take the model Hamilto- 
nian of the system as iJ = ^ Hj, where j is the layer 
index. 
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where Vho = l/2Afw^r2 + l/2Muj'^^a'^f is the external 
trapping potential with r^ = x^ +y'^, V2D is the 2D gra- 
dient operator and U2D is the 2D interaction strength. 
The operator tp]^{r) creates a fermion of mass M in in 
jth plane with spin a =t,i at position r — {x,y). We 
consider a tight ID lattice where interlayer tunneling en- 




the chemical potential. One can treat this local chemical 
potential with local density approximation. In this pa- 
per, we take the external trapping potential into account 
when we calculate the dipole oscillations, otherwise we 
treat the system as homogenous. In order to take into 
account superfluid flow along z-direction, we take super- 
fluid order parameter in the form A,„ — A exp[iTO0] . The 
periodicity along z-direction allows us to write the Fermi 
operators. 



FIG. 1: Tunneling energies as function of laser intensity. Gray 
line: Tunneling energy (tt) proposed in ref '§| for 3D optical 
potentials. Black line: Tunneling energy (ta) proposed in 
ref. [3] for ID optical potentials. The dimensionless parameter 
s accounts the laser intensity. 



ergy t is small and consider the following expression pro- 
posed in ref. [7| and ref. [8| for the weakly coupled layers. 
In terms of dimensionless parameters, the tunneling en- 
ergy for a ID lattice proposed in ref. [7| is ta/Ej^ = 
7Ts{Tr'^ /A — 1) exp[— 7r^^/7rs/4]. This contrasts with the 
isotropic three-dimensional (3D) lattice tunneling param- 
eter tb/En = (2s3/4/0f) exp[-2Vs] given in ref. [l]. 
Both of these are derived using a harmonic approxima- 
tion around the minima of the optical lattice potential. 
Notice that these ratios can be varied by changing the 
laser intensity. As depth of the optical potential in- 
creases, the atomic wave function becomes more and 
more localized and the tunneling amplitude decreases. 
For comparison, these tunneling amplitudes are shown 
in FIG. ([1]) as a function of the laser intensity. For 
asymptotically deep lattices and any separable poten- 
tial, the tunneling energy should be independent of the 
dimensionality of the lattice. The asymptotically differ- 
ence between ta and tb seen in FIG. (P) may be due to 
the neglecting of the non-orthogonality of the adjacent 
lattice Wannier functions in ref. [71]. In this paper, we 
use the expression t = tb for all our calculation. In the 
limit t ^ 0, the system is decoupled planes of Fermi 
atoms. We use mean field theory to decouple the in- 
teraction term writing f^2D'0lt(f")V'Ii ('')V'ji('')V'jT('') ~ 
[A](r)VjiWV'jTW + h.c] - \Aj{r)\yU2D- The Fourier 
transform of the Hamiltonian gives. 
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where e^ = n^k'^/{2M) with k"^ = kl + k"^. Notice that 
we have absorbed the external trapping potential into 
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Transforming the Hamiltonian given in Eq. ^ using 
above transformation followed by the usual Bogoliubov 
transformation, the Hamiltonian per plane can be ex- 
pressed as |9| 
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where N is the number of planes and Ck — Cfe + 
2tcos(A:2d)cos((/)/2). Notice that the quasi particle en- 
ergy depends on the phase of the superfluid order pa- 
rameter. In the normal phase, the phase factor (jj can 
have multiple values of tt [10|. The grand potential of 
the system O = (-1//3) ln[ZG] with Zq = Trcxp[-f3H] 
at zero temperature is 
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The 2D contact interaction U2D is related to the bound 
state energy Eb as [13| 
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The bound state state energy Eb = 

(CTiWi/Tr) exp[\/27r/L/as], where Os is the 3D s- 
wave scattering length and C w 0.915. Notice that the 
bound state energy depends not only on the scattering 
length but also on the lattice potential. Performing the 
momentum integrals, the grand potential is 
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FIG. 2: Superfluid order parameter A as a function of Ea ~ 
h^ /{2Mal) for as < (long dashed line), as -^ oo (short 
dashed line) and as > (solid line) at fixed density n — 
2.5m/{2-Kh'^) and s = 8. 



Josephson current and Bloch oscillations 

For the solid-state electronic systems, the size of the 
Cooper pair is usually larger than the lattice spacing and 
one can treat the system as homogenous when calculat- 
ing physical quantities such as the supercurrent. How- 
ever, the periodicity of the cold atomic system is ad- 
justable. As a result one can expect new phenomena 
on optical lattices with Fermi superfluid. For example, 
oscillating fermionic superfluidity which is absent in the 
bulk solid can be observed with the periodic potentials. 
These oscfllations may be used as a probe of BCS-BEC 
crossover and of the normal to superfluid transitions. We 
calculate the particle supercurrent along z-direction us- 
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Then the gap equation and the number density are cal- 
culated by using the equations, dfl/dA — and n — 
—dfl/dfi, respectively. 
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III. RESULTS AND DISCUSSION 

All our derivations for physical quantities will be cal- 
culated from Eq. dZ]), Eq. dH) and Eq. ^. For example, 
the solutions of Eq. ^ for various parameters are shown 
in FIG. ^. For the case of decoupled layers (in the limit 
t = 0), Eq. dl]) gives A == ^/Eb{Eb + 2^i). This results 
the grand potential Q = —[m/{8nti^)]{EB + '2ii)^ and the 
number density n = ni/{2Trfi )]{Eb + ^ji). 

In the following subsections we calculate the Joseph- 
son current, the sound velocity and the dipole oscilla- 
tion frequencies. First we calculate the oscillatory super- 
fluid velocity along the direction of optical lattice through 
the Josephson current in the absence of harmonic poten- 
tial. Next, we consider the propagation of sound along 
z-direction and calculate the sound velocity in a homoge- 
nous system. Later we include the external harmonic po- 
tential to investigate the dipole oscillations of the atomic 
cloud. 
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In early quantum theory of electrical conductivity, 
Bloch and Zener [ll| predicted that when a static electric 
field is applied to a crystalline electronic system, instead 
of uniform motion one would naively expect, electron 
in the crystal should oscillate. However these so-called 
Bloch oscillations have never been observed in the bulk 
crystals. The reason is that the period of Bloch oscil- 
lations is much larger than the electron scattering rate. 
As the cold atomic superfiuid in optical lattices are clean 
and controllable, these systems have shown to be well 
suited for the observation of Bloch oscillations. In fact, 
the oscillations are already observed in one dimensional 
optical potentials |12l |. 

We use superfluid velocity along the z-direction to get 
some understanding about the amplitude of the Bloch os- 
cillations. Consider an atom is in the Bloch state \l,kQ >, 
where I is the discrete band index and fco is the continu- 
ous quasi momentum. If the atom feels a constant force 
F (In optical lattices, one can mimic a constant force on 
atoms by introducing a tunable frequency difference be- 
tween two counter propagating laser beams), the above 
Bloch state evolve to a another state \l,k{T) >, where 
fc(T) = ko + Fr/h is the quasi momentum at time r. 
We assumed that the force is weak enough and the ID 
lattice is tight enough not to induce the inter band tran- 
sitions. Therefore setting cf> = Ft/ ft in Eq. pTj) . we 
can calculate current due to the induced Bloch oscilla- 
tions. Then we calculate the superfluid velocity using 
vsF[FT/h)=3{FT/n)/n. 
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FIG. 3: Block oscillations in (a). BCS regime, (b). Unitar- 
ity regime, and (c). EEC regime, where we fixed the scat- 
tering length to be a^ — —^/h^ /{2mEii), as -^ oo, and 

Us — -\-\Jf^ l{\lmEn) respectively. The parameter s varies as 
4 (black), 5 (dark gray) , 6 (light gray) and, 8 (lighter gray). 
For aU the calculation, we use the density n — 2.bm/{2-K'h ). 



A'^cos'^ {Ft /2h)t^ + 0{t'^) 



These superfluid velocities are shown in FIG. ([3]) for 
three different regimes. It is clear that the Bloch oscil- 
lations are stronger (the amplitude is larger) in the BCS 
regime where A is smaller and where the coherence length 
is large compared to the inter layer distance. In other 
words, away from the BCS regime the effective mass of 
the Cooper pair becomes larger (due to the strong inter- 
action) and localizes in each plane. In order to observe 
these Bloch oscillations, superfluid velocity should not be 
exceeded the critical velocity. This Landau criterion for 
superfluidity is satisfied if the Cooper pair size is smaller 
than the lattice periodicity. 



Propagation of sound 

We consider a homogenous system and the prop- 
agation of the sound along z-direction. The form 
of the phonon dispersion lo = CzQ gives the sound 




FIG. 4: Sound velocity along the direction of ID lattice as 
a function of dimensionless parameter s. For small tunneling 
energies (tight ID potentials), the sound velocity is almost 
the same in the entire BCS-BEC crossover region. 



velocity along z-direction c^ = y {n/M){dfi/dn), 

where the effective mass of the atoms M is given 
by Af-i = (d^Et/dk^), where Et = 2t[l - 
cos{k,d/h)]. Using (d^Et/dk^) = 2t(P /n^{cos{k,d/K)) = 
{2tcP/h^{{l/n)^f^ f, nkCos{kzd/h) with atom density 
n ~ ^j, J, Uk and summing over the momentum states, 
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Notice that the presence of periodic lattice is taken into 
account through both the effective mass and the com- 
pressibility dfi/dn. And the ratio c^/cl = M, where c± 
is the sound velocity in transverse direction. Calculating 
the compressibility using the number equation (Eq. ([§]), 
the sound velocity is 
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The sound velocity as a function of s for = is given 
in FIG. ([1]). As seen in Eq. (fT^]) . the sound velocity is 
independent of A and /x for the lowest order in tunnel- 
ing t. As a result, the sound velocity is almost same in 
the entire BCS-BEC regime for a very tight ID optical 
lattice. In fact, this result holds for a non- interacting 
degenerate Fermi gas in a very tight periodic potential. 
In the absence of a periodic potential, the sound velocity 
is smaller than the Fermi velocity. Therefore the sound 
cannot propagate in a degenerate non-interacting Fermi 




m_i = No/{2Miu'^) can be calculated. Here H is the 
Hamiltonian, [A, B] = AB — BA is the usual commuta- 
tor, Nq is the total number of atoms, and (..) is average 
over the ground states. Then the upper bound of the 
dipole oscillation frequency lu^ is calculated from the ex- 
pression 
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This gives the frequency of the center of mass oscillation 

along the ID lattice uJd = uJzy M/M. Notice in the 
absence of optical lattice LOd = i^z is independent of the 
interaction. 

The dipole oscillation frequency as a function of s is 
given in FIG. ^ for = 0. In the entire BCS-BEC 
region, the dipole oscillations show significant reduction 
in the presence of ID lattice. This reduction is stronger in 
the BEG regime where A is larger. As evidence from the 
FIG. (O, the oscillation frequency shows a weak chemical 
potential (hence the total number) dependence. 



IV. CONCLUSIONS 



FIG. 5: Dipole oscillations in (a). BCS regime, (b). Unitar- 
ity regime, and (c). BEC regime, where we fixed the scat- 
tering length to be a^ — — y ?i^/(2m_E_R), as -^ oo, and 

fls = +yj%^ /{2mER) respectively. The parameter ^ varies 
as 0.5-Bfl (black), 0.75-Efl (dark gray) and. En (light gray). 



gas. However, the periodic potential modifies the effec- 
tive mass of the atoms and the sound can propagate even 
in a non-interacting Fermi gas. 



We study some of the dynamical properties of a two- 
component Fermi gas trapped in a parabolic plus ID pe- 
riodic potential. We use a zero temperature mean field 
theory in the limit of tight ID confinement. We find that 
the Bloch oscillations become weaker as the interaction 
strength changes from weakly attractive limit to weakly 
repulsive limit through the unitarity regime. Further- 
more, we find that the sound velocity is independent of 
the interaction for very tight ID confinements. Finally, 
we find that the dipole oscillation frequencies strongly 
depend on both the interaction and the ID periodic po- 
tential. 



Dipole oscillations 

The Bloch Oscillations discussed in the previous sub- 
section are due to the atoms oscillations along the Bloch 
band under a linear force, and the dipole oscillations or 
the center of mass oscillations are due to the confinement 
in the external trapping potential. Even though the the- 
oretical difference is clear, the probing of these two types 
of oscillations may be very complicated. One possible 
way of distinguishing these two types of oscillations is 
from observation of the damping rate. If the momentum 
distribution is narrower than the first Briloin zone, then 
the damping rate of the Bloch oscillation will be zero. 

The frequency of the dipole oscillations can be derived 
from the sum rule approach. Defining the dipole op- 
erator D = ^ Zj, the energy weight moment mi = 
(1/2) ([_D, [H, D]]) and the inverse-energy weight moment 
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